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Spin qubits composed of either one or three electrons are realized in a quantum dot formed at a Si/Si02 
interface in isotopically enriched silicon. Using pulsed electron spin resonance, we perform coherent control 
of both types of qubits, addressing them via an electric field dependent ^-factor. We perform randomized 
benchmarking and find that both qubits can be operated with high fidelity. Surprisingly, we find that the g- 
factors of the one-electron and three-electron qubits have an approximately linear but opposite dependence as a 
function of the applied dc electric field. We develop a theory to explain this ^-factor behavior based on the spin- 
valley coupling that results from the sharp interface. The outer “shell” electron in the three-electron qubit exists 
in the higher of the two available conduction-band valley states, in contrast with the one-electron case, where 
the electron is in the lower valley. We formulate a modified effective mass theory and propose that inter-valley 
spin-fiip tunneling dominates over intra-valley spin-fiips in this system, leading to a direct correlation between 
the spin-orbit coupling parameters and the ^-factors in the two valleys. In addition to offering all-electrical 
tuning for single-qubit gates, the ^-factor physics revealed here for one-electron and three-electron qubits offers 
potential opportunities for new qubit control approaches. 


Silicon is known among the semiconductors to have small 
spin-orbit coupling (SOC), a beneficial fact for silicon quan¬ 
tum computing, since charge noise is largely decoupled from 
information stored in the spin m. Furthermore, silicon can be 
isotopically enriched and chemically purified to ^^Si, thereby 
removing nuclear spin background fluctuations and so silicon 
is often referred to as a semiconductor vacuum ||2|- These two 
facts have motivated intense research on silicon qubits, lead¬ 
ing to recent realizations of single-qubit 130 and two-qubit 
logic gates. Despite the small SOC in isotopically purifled 
silicon quantum dots, the small tunability of the ^-factor via 
gate-controlled electric flelds allows one to electrostatically 
turn on and off the spin rotations that constitute single-qubit 
gates (30, thereby providing an important tool for quantum 
computation. 

The low-energy subspace in silicon quantum dot (QD) sys¬ 
tems is governed by two spin-degenerate valley states. When 
these valley states are quasi-degenerate, qubit operation be¬ 
comes complex (6), and the coupling of qubits is even more 
challenging cni. However, the valley states can be separated 
using a vertical electric fleld and the sharp potential of an in¬ 
terface (3 El, and their energy separation can be electrically 
controlled over several hundreds of /ieV. While one-electron 
spin qubits are naturally operated in the lowest valley state 
(313, it is intriguing to consider the potential performance of 
qubits operated in the higher valley state, which is also known 
to have a long spin lifetime if orbital relaxation can be sup¬ 
pressed m. When spatial conflnement in the QD is strong, 
the orbital excited states are lifted high in energy and qubit 
operation in the upper valley state is possible by populating 
three electrons in the quantum dot. In this mode two electrons 
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FIG. 1: (a) Schematic representation of the quantum dot sys¬ 

tem. The quantum dot is defined using the confinement gate CG 
and plunger gate PG and the yellow shading represents the regions 
where electrons are accumulated, with Fz the perpendicular electric 
field direction. ESR control is via a dc magnetic field Bq = 1.4 T 
(in the Si [110] in-plane direction) and an ac magnetic field Bi. (b) 
The quantum dot qubit can be operated using the spin states of one 
electron, or using three electrons, where two electrons (blue) occupy 
the lowest energy valley state and the third electron (red) is in the 
higher energy valley state, (c) Charge stability diagram showing the 
electron occupancy N in the quantum dot, measured with a nearby 
SET. 


form a singlet in the lower valley state and the third electron 
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is operated in the upper valley (see Fig[^ and[^). It has been 
suggested that such multi-electron qubits could enhance the 
gate fidelity, due to partial screening of electrical noise ica. 


In this letter, we demonstrate high-fidelity operation of one- 
and three-electron spin qubits, operated in the lower and up¬ 
per valley, respectively. Using pulsed electron-spin-resonance 
(ESR) we map out with high precision the qubit frequency as a 
function of the applied perpendicular electric field, controlled 
with electrostatic gates. We experimentally demonstrate and 
theoretically explain how inter-valley spin-orbit coupling at 
the Si/Si 02 interface results in an opposite dependence of the 
^-factor for the two valleys, that is correlated with the cor¬ 
responding dependence of the 2D Rashba and Dresselhaus 
spin-orbit coupling in each valley. We also present random¬ 
ized benchmarking on both qubit systems, showing they are 
capable of fidelities above 99%, approaching the surface code 
thresholds for fault-tolerant quantum computing ca. 

The QD structure is fabricated on an epitaxially grown, iso- 
topically purified ^®Si epilayer with a residual concentration 
of 29si at 800 ppm E). The aluminum gates are fabricated 
with electron-beam lithography using a multi-level gate stack 
silicon metal-oxide-semiconductor (SiMOS) technology ifT^ . 
see Fig. The charge stability diagram of the quantum dot 
is shown in Fig. [^. We use ESR to control the one-electron 
qubit El and the three-electron qubit, see Fig. From a 
Ramsey sequence we find a long dephasing time T 2 = 70 /is, 
which is slightly less than we have previously measured for 
the one-electron qubit, which had T 2 = 120 /is El 

We have demonstrated electric field control over the res¬ 
onance frequency z^esr of the one-electron qubit f7|, show¬ 
ing tunability over several MHz that appears linear in elec¬ 
tric field, corresponding to more than 3000 times the 2.4 kHz 
ESR line-width. We find that spin-valley mixing of the 
QD eigenstates due to interface roughness (TTIl would pre¬ 
dict a modification of the electron ^-factor that is two or¬ 
ders of magnitude smaller than is found experimentally, to¬ 
gether with a non-linear dependence close to the anticross¬ 
ing point of the spin-valley states that we do not observe. 
Here we propose and analyze a model where the ^-factor 
modification proceeds via inter-valley spin-flip tunneling, me¬ 
diated by the strong ^-confinement at the interface. The 
Si/Si02 (001) interface of silicon MOS quantum dots can 
be described with a Hamiltonian that consists of a bulk 
term Hq and an interface term Hif. The reduction of the 
bulk Si crystal symmetry at the interface, in the presence of 
strong perpendicular electron confinement induced by an ap¬ 
plied electric field lifts the six-fold valley degeneracy, 
leaving two low-lying A-valleys at Fkoz. These are then 
mixed via enhanced inter-valley tunneling due to the strong 
^^-confinement at the interface CSIISl. The consequent effec¬ 
tive two-valley Hamiltonian acts on the four-component vec- 
tor where 
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FIG. 2: Demonstration of qubit control of the three-electron qubit, 
(a) 2D colour map showing Rabi control of the spin-up fraction /^, 
by varying the microwave pulse length and the microwave driving 
frequency vesr- We have subtracted a reference frequency rq — 
39.045 GHz (corresponding to ^ = 1.9908) for clarity, (b) Ramsey 
fringes, obtained by varying the waiting time in between two ESR 
7r/2 pulses. The decay in the spin-up fraction corresponds to T 2 
= 70 ps. The confinement gate voltage is Vc = -0.2 V. 


the bulk part (spin and valley degenerate) is given by 
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with the quasi-momentum operators kj = —idj', and Ux,y = 
+ y‘^) and Uz = \e\FzZ are the in-plane and per¬ 
pendicular confinement electron potentials, respectively. Here 
mi, mt are the Si effective A-valley electron masses, |e| is the 
electron charge, and fi is the reduced Planck constant. Taking 
into account the large band offset of Si/Si02, the interface 
term is 


~ ^ 0 ) - ii :,—- zo)kz 

2Rmi 2mi 

+ S(z - zo)Vi{(k), 


( 2 ) 


where is a parameter with dimension of length, character¬ 
izing an abrupt interface |[l7][l8l, and \R\ Iz <C Id’, here 
Iz = /2mi\e\Fz)^^^ and Id = (fi/mtcco)^/^ are the per¬ 

pendicular and in-plane confinement lengths (assuming much 
stronger ^-confinement). For R ^ 0 the interface Hamilto¬ 
nian (Eq. corresponds to the standard infinite boundary 
condition (BC) ^{z) lz=zo= 0^ while for finite R it generates 
spin and valley mixing at the interface, z ^ zq. Following the 
symmetry reasoning of Refs. HMD the spin-valley mixing 
interface matrix Vif(k) can be expressed via the C 2 V invariants 
Hji(k) = a^zky - cFykz;, Hoik) = a^k^z - o-yky, resulting in 


Vidk) 


A{k) Vh + B{k) 
V*/2 + -Bt(fc) A{k) 


( 3 ) 


In Eq. © the 2x2 block-diagonal element A{k) = 
SDHD{k) - 1 - SRHR{k) corresponds to intra-valley spin- 
flipping transitions, while the off-diagonal elements VI 2 , 
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B{k) = XDHD{k) + XRHR{k) are related to inter-valley 
tunneling (in momentum space) with no spin-flipping or with 
a spin-flipping process. Since experimentally the valley split¬ 
ting ( r\j I y |) is generally large with respect to the spin-flipping 
terms, it makes sense to diagonalize with respect to the lead¬ 
ing V -matrix element and via a unitary transformation we find 


vy{k) = 


\V\ + A+\Ba \B,s 

h.c. -\V\+A-\Ba 


(4) 


This matrix is approximately diagonal in the valley basis 
l^i), 1 ^ 2)5 with a calculated valley splitting energy T^vs = 
2\V\R^\ip'{{))\^ = 2\V\R^l~^ (X Fz. We neglect the off- 
diagonal contribution = B — in Eq.(|^, since it 

is suppressed as ^ 1/E’vs and E’vs is typically several hun¬ 
dreds of /ieV in MOS quantum dots □ [m. Thus, in the 
valley subspaces \vi), \v 2 ), one can consider two independent 
boundary conditions as in Eq.Q, with spin-flipping interface 
matrices Ki,v 2 = ^ T ^ T \ 

in which the inter-valley spin-flip tunneling element changes 
sign between vl and v2. 

The effective 2D spin-orbit Hamiltonians (proportional to 
the Rashba and Dresselhaus forms, HR{k), HD{k)) are 
straightforwardly calculated from Eq.Q, with the corre¬ 
sponding 2D SOC parameters changing sign as well: 

C^R-^vl,v2 = 2[Si? ^ \xr\ COs(0i? - 0y)] R^\ip'{0)\‘^ 

Pd-v1,v2 = ‘^[sd T \Xd\ COs{(I)d — (pv)] (5) 


The scaling of the spin-orbit terms with the electric field 
is linear, as is the valley splitting, E^vs- Here, we have intro¬ 
duced the phases (pR and (j )d for the Rashba and Dresselhaus 
terms and (/:?'(0) is the derivative of the 2 ;-component of an 
eigenstate of the bulk Hamiltonian Rq. These results are sim¬ 
ilar to the strong field limit results of Ref. ED 

Explicit calculation of the ^-factor change, based on the 
interface Hamiltonian Eq.Q and the fact that an in-plane 
magnetic field mixes the perpendicular and in-plane motion, 
shows that the in-plane ^-factor renormalization dg\^ propor¬ 
tional to aR, Pd. Eor the magnetic field parallel to the [110] 
direction (as in the experiment) one obtains the expression: 

r {OLR.yly2—PD-vl,v2)\e\, \ 

^gvXv 2 = --r-- {z) (6) 

hllB 

where /ie is the Bohr magneton and (z) 1.5587/^ is an 

average of the 2 ;-motion in the lowest subband, see Eq.Q. 
The ^f-factor scales as Fz^^, which is close to a linear scaling 
over the range (^ 10%) of the experimentally applied electric 
fields, see Eig. 3b. 

We therefore expect from Eq.([^ that the renormalization 
6g will be of opposing sign for the two valleys, following the 
sign change of the SOC parameters in Eq.Q. In particular, the 
change will be exactly opposite for zero intra-valley spin-flip 
coupling, Sr, sd = 


(7) 


b |EJ (MVm-^) 

27.5 28.0 28.5 29.0 29.5 



FIG. 3: (a) Magnetic field dependence of the resonance frequency 
of the one- and three-electron qubit, with Fz = 28.25 MV/m. The 
experimental data of both qubit systems has been subtracted by the 
^le-factor for comparison and we have calibrated the dc magnetic 
field using the crossing point of the one- and three-electron qubit res¬ 
onance frequencies. We find gie = 1.9975 and gse = 1.9912. (b) Gate 
tuned electric field control over the valley ^-factor at Bq = 1.4015T. 


Relatively smaller corrections due to non-zero intra-valley 
spin flipping, sr,sd % 0 will generally violate Eq. Q, leav¬ 
ing the ^-factor changes opposite in sign, but with different 
absolute value, \6gyi\ % \Sgv 2 \- 

To observe this experimentally, we control the quantum 
dot electric field via the plunger gate PG and the confine¬ 
ment gate CG, see Eig. In Eig. we show the mag¬ 
netic field dependence and in Eig. we show electrical 
control over the qubit resonance frequency vesr^ The op¬ 
posite electric field dependence of the ^-factor for the two 
valleys is in qualitative agreement with the prediction of 
Eq.0. Since the resonance frequency of the one-electron 
qubit increases with the electric field, while the resonance fre¬ 
quency of the three-electron qubit decreases (see Eig. |3>), 
we infer from Eq.([^ that the Rashba and Dresselhaus con¬ 
tributions are in this experiment subject to the constraints: 

<^Xinter-val = \Xr\ COs(0i? - py) “ \Xd \ COs{(I)d “ Pv) > 0, 

^^5intra-vai = Sr — Sr <0. The change in sign of 6g is evi¬ 
dence that the inter-valley spin-flip contributions dominate the 
intra-valley spin-flip processes and from the J^-dependence 
we estimate the ratio ^i^Xinter—vai/l^'^intra—vai| ^ 2.6. This 
observation is consistent with tight-binding calculations on 
SiGe quantum wells ED, which predict that the inter-valley 
transitions can be about an order of magnitude larger than the 
intra-valley transitions. Interestingly, the values and signs of 
the SiGe parameters, as substituted in Eqs. 0 and ([^, repro¬ 
duce also the correct qualitative behavior of (Sgy 2 ) that 
increases (decreases) with the applied electric field Fz, as is 
observed experimentally in Eig. [^. However, the experimen¬ 
tal ratio of the ^f-factor changes is \Sgy 2 \/\Sgyi\ ^ 2.2, while 
that calculated with the Si/SiGe parameters is ^ 1. Such dif¬ 
ferences can be expected due to the much greater band-edge 
offset in Si/Si02, disorder 1^ . and built-in electric fields 
which may also influence the theoretical results. 

In order to explore the qubit performance for quantum com- 


Sgvi = -Sgv2- 











4 



b RB sequence ^ 



Random 

Clifford 

M Reeovery 
l^ln Clifford 

c 

Qubit fidelity 


3e4 

F. 

Referenee 

99.0% 

97.6% 

Average H 

99.1% 

97.2% 

F, . 

Referenee 

99.9% 

99.3% 

Average H 

99.9% 

99.3% 


FIG. 4: Clifford based randomized benchmarking (a) Sequence fi¬ 
delity as a function of the sequence length and (b) schematic repre¬ 
sentation of randomized benchmarking, where H is an interleaved 
test gate. In (a), the black filled circles correspond to standard ran¬ 
domized benchmarking and the green open circles to the average 
of the single qubit gates [/, zbX, =b|X, zbY, =t|y], obtained by in¬ 
terleaved randomized benchmarking and normalizing the sequence 
length with 2.875/1.875, such that it matches the average Clifford 
gate length. Both data sets are fitted with a two-fidelity model (see 
text) and the results are shown in (c), where the standard error is 
smaller than the corresponding gate error. 


putation, we have performed (interleaved) randomized bench¬ 
marking (RB) |[23l[24]| on the one-electron qubit 0 and three- 
electron qubit, and all results are shown in Fig. In order 
to eliminate the free fitting parameter B, which is a constant 
offset parameter present in standard RB fits, we plot the se¬ 
quence fidelity combination F = — 1, which ap¬ 

proaches zero for infinite sequence length when the assump¬ 
tions of RB hold ll^ . When the noise is gate independent, 
an exponential decay is expected for both standard and in¬ 
terleaved randomized benchmarking. However, when low- 
frequency noise is present, non-exponential decays arise 1 ^ . 
This non-exponential decay is due to slow drifts in the reso¬ 
nance frequency, such that the time ensemble is averaged over 
sequences with small detuning (resulting in a high fidelity, 
Fhigh, and a slow exponential decay) and large detuning (re¬ 
sulting in a low fidelity. Flow, and a fast exponential decay). 

When such low-frequency noise is present, the fidelity 
varies over time and we use a two-fidelity model to analyse 
the data 1^ . We have fitted the data using F = 
where A quantifies the SPAM error and p and q are two polar¬ 
ization parameters, and in Fig. we show the corresponding 
fidelities. The three-electron qubit has a relatively low fidelity 
when the noise causes a large detuning {Flow ~ 97%). How¬ 
ever, when the microwave driving frequency is on resonance, 
both qubits have a fidelity above certain thresholds for fault 
tolerant quantum computing lEl; the average single-gate fi¬ 
delity being F^igh = 99.9% for the one-electron qubit system 
and Fhigh = 99.3% for the three-electron qubit system. While 
the three-electron qubit initially shows a non-exponential de¬ 
cay, for higher m the decay approaches a pure exponential, 
indicating that low frequency noise has little impact in this 
range. Since both qubit systems are operated with the same 
setup, we expect similar calibration errors. This exponential 


decay is therefore likely due to high-frequency noise. 

The faster decay of the sequence fidelity of the three- 
electron vs. one-electron qubit is consistent with a larger sen¬ 
sitivity to electrical noise, as revealed by the larger frequency 
shift with gate voltage, \6i'^e/^V\ 2.2\di'ie/6V\, shown in 
Figj^. The frequency detuning caused by electrical noise re¬ 
sults in rotations around the 2 ;-axis of the qubit Bloch sphere 
and opposite to the Rabi driving axis. By taking the small- 
angle approximation, we find that the gate error increases with 
the square of the noise term. This results in an error rate that is 
around 5 times larger for the three-electron qubit, comparable 
with the difference in fidelities between the one-electron and 
three-electron qubits. It is therefore likely that both qubits are 
ultimately limited by high-frequency electrical noise. 

The recent realizations of single- and two-qubit gates using 
isotopically purified silicon quantum dots EEl are now re¬ 
vealing the early promises of silicon as a platform for quantum 
computation and the possibility of qubit operation with either 
one-electron or three-electrons allows more flexibility in scal¬ 
ing these systems. The ultra-narrow spin-resonance linewidth 
associated with the long coherence times of these qubits has 
pushed silicon into a new regime, where the weak spin-orbit 
coupling in silicon becomes not only visible, but also forms 
a new tool to control the spin states, as shown here. The re¬ 
markably large electric fleld control in SiMOS quantum dots 
provides further motivation to explore spin-orbit coupling in 
silicon for qubit control and spin manipulation. 
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